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Abstract 

Fractional (in time and in space) evolution equations defined on Dirichlet regular bounded 
open domains, driven by fractional integrated in time Gaussian spatiotemporal white noise, 
are considered here. Sufficient conditions for the definition of a weak-sense Gaussian solution, 
in the mean-square sense, are derived. The temporal, spatial and spatiotemporal Holder 
continuity, in the mean-square sense, of the formulated solution is obtained, under suitable 
conditions, from the asymptotic properties of the Mittag-Leffier function, and the asymptotic 
order of the eigenvalues of a fractional polynomial of the Dirichlet negative Laplacian operator 
on such bounded open domains. 
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1 Introduction 

Space-time fractional diffusion equations are introduced when integer-order derivatives in space 
and in time are replaced by their fractional counterpart (see, for example, [43] ). In particular, they 
can model anomalous diffusion processes in physics (Meerschaert et al. [35]). Fractional diffusion 
equations are very popular in several fields of application (see Gorenflo and Mainardi [53] ; Metzler 
and Klafter [41], among others). 

The development of space-time fractional diffusion equations is supported by the theory of 
operators for generalized fractional calculus (see, for example, Kiryakova [53|, [50]; Machado, 
Kiryakova and Mainardi [3S] , [31] , and the references therein). 

Since the pioneer papers by Bochner m and Feller m who proved the connection between 
the stable distribution and fractional calculus, the theory of a-stable distributions and processes 
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has been extensively developed. Specifically, Bochner m formulated the Cauchy problem, whose 
solution is the symmetric a-stable distribution. Feller [20] extended these results to a more general 
situation by replacing the fractional Laplacian — (—by a pseudodifferential operator with 
symbol 

— |A|“ exp (i sign(A) 07r/2), A G R, a £(0,2), 

where a is the index of stability, and 9 is the index of skewness (asymmetry). The corresponding 
solutions generate all stable distributions. The Cauchy problem characterizing the main properties 
of the Riesz-Bessel distribution, generated by the pseudodifferential operator (—A)“/^(J — A)'>'/^, 
was analyzed in Anh, Leonenko and Sikorskii |^. In spite of a large number of fractional oper¬ 
ators (see Samko et al. |45] 1. there were few known specific examples producing more general 
distributions. 

The traditional model for particle spreading at the macroscopic level is the well-known heat 
equation 

dtu = Au, 

with A denoting the Laplacian operator and dt the partial derivative in time. The relative particle 
concentration can be predicted in terms of the Gaussian probability density providing a point 
source solution of the heat equation. The paths of individual particles are described in terms of 
the realizations of Browninan motion. Dirichlet boundary value problems for the heat equation, 
as well as for more general equations, given in terms of elliptic diffusion operators, are detailed 
in Bass m and Davies among others. The phenomena of particle sticking and trapping can 
be described by using the fractional derivative df for 0 < /3 < 1 instead of the partial derivative 
in time dt- On the other hand, long particle jumps can be represented by the fractional power 
(—A)“/^, for 0 < a < 2, in place of the negative Laplacian (—A). The space-time fractional 
diffusion equation is then defined in terms of both fractional derivatives in time and in space: 

af u = (-A)“/ 2 u, ( 1 ) 

whose solution displays self-similarity and heavy tails. The particle concentration profile pro¬ 
vided by the corresponding probability density solution has sharper peak and heavy tails. A 
non-Markovian setting can then be introduced through time change by an inverse stable subor- 
dinator (see also Barkai et al. [5]; Benson et al. [a; Gorenflo and Mainardi jUj; [53]; Meer- 
schaert et al. ESI; Schneider and Wyss [33], among others). An extension to the case of Riesz- 
Bessel subordinators was addressed in Anh and McVinish [7] , using the pseudodifferential operator 
(—A)“/^(/—A)'’'/^. A different stochastic framework is analyzed in the papers by Anh and Leonenko 
m, 0 , where the spectral representation of the mean-square solution of the following stochastic 
space-time fractional equation with random initial conditions 

5fu = (-A)“/2(j_ a)^/2u, uo(x)=r?(x), x £ R", 

is derived. Here, 77 is a measurable random field defined on a complete probability space (D, A, P). 
Gaussian and non-Gaussian limiting distributions of the renormalized solution are obtained as 
well. Also, in the context of stochastic evolution equations on an unbounded domain, a functional 
approach was adopted by Kelbert, Leonenko and Ruiz-Medina EH], where the spectral properties 
of the mean-square solution of fractional in time and in space evolution equations driven by random 
white noise are derived. These results are extended to the more general framework of stochastic 
partial differential equations driven by fractional Brownian motion in Leonenko, Ruiz-Medina, 
Taqqu m, where the spectral correlation structure of the mean-square solution to fractional space- 
time evolution equations driven by fractional Brownian are analyzed. In the particular case of the 
domain being an n-dimensional rectangle, starting from the pioneering papers on the stochastic 
heat equation by Walsh [48], [49], Angulo, Ruiz-Medina, Anh and Grecksch [3] considered the 
following fractional (in space) version of the heat equation: 

(t, x) + (-A)“/^ (/ - A)'^/^ c (t, x) =e {t, x) , ( 2 ) 
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for t G K+, X G D = {0, Li) x • • • x (0, Ln), and n<a + j<n + 2, where e {t, x) is a zero-mean 
Gaussian space-time white noise. The space-time mean-quadratic and sample-path local variation 
properties of the solution are derived in this paper. The unbounded domain case is addressed as 
well. Angulo, Anh, McVinish and Ruiz-Medina [5] derive the extension of these results to the 
case of fractional derivatives in time and space with the same spatial pseudodifferential operator 
(—(/ — and D = (0,Li) x • • • x (0,L„) (see also the references therein). 

In the context of fractional diffusion on bounded domains, we refer to the papers by Defterli, 
D’Elia, Du, Gunzburger, Lehoucq and Meerschaert [18]; Chen, Meerschaert and Nane [M], and 
Meerschaert, Nane and Vellaisamy [31], where strong solutions, and their probabilistic representa¬ 
tion are obtained. On the other hand, Mijena and Nane [H] consider fractional heat equation on 
unbounded domains, with a non-linear random external force, involving space-time white noise. 
Sufficient conditions for the existence and uniqueness of mild solutions, as well as for their conti¬ 
nuity are derived. We consider here a different framework. Specifically, we study the weak-sense 
solution of the following fractional in space and in time stochastic partial differential equation, 
with Dirichlet boundary conditions, and null initial condition; 

-^c {t, x) -k (/ - A^,)'*'/^ c {t, x) = ll~^e (t, x), kgD (3) 

c(<,x)=0, xGdD, Vt, c(0,x) = 0, Vx G D C K", (4) 

for /I G (0,1), a -|- 7 > n, where equality is understood in the mean-square sense. Here, the driven 
process 

u)-f^e{u)du (5) 

is constructed from space-time Gaussian white noise e, defined on a basic probability space (D, A, P), 
and satisfying 

£;[e(t,x)e(s,y)] = d{t - s)S{x - y), 

for all t,s G M+, and x, y G D, with S being the Dirac Delta distribution. Specifically, the driven 
process is constructed from fractional integration of order /3 — 1, in time, of the space-time Gaussian 
white noise e, where integration is understood in the mean-square sense (see, for example, Samko 
et al. m)- It is well-known that the inverse of the fractional integral of order /3 — I is the fractional 
derivative of order 1 — /3. This fractional derivative, considered in this paper, is the regularized 
fractional derivative in time or fractional-in-time derivative in the Caputo-Djrbashian sense: For 
/3e(0,i], 


dtP I (t-r) '^u(T,x)dr-if/3 G (0,1), t G (0,T] 

(see Meerschaert and Sikorskii [40]; Podlubny [44]). 

Although we refer here to the particular case where fractional derivatives in space are defined 
from the operator (—Ad)“^^ (/ — A^)^^^ , with (—An) representing the Dirichlet negative Lapla- 
cian operator on regular bounded open domain D, the results derived in this paper hold, in general, 
for a fractional polynomial of the Dirichlet negative Laplacian operator on D, with constant coef¬ 
ficients, as proved in Theorem 18.11 in Section [5] In this paper, special attention has been paid to 
operator (—Ad)^^^ (J — A/j)^^^ , since, for suitable domains, e.g., for bounded open domain satis¬ 
fying the exterior cone condition, the eigenvalues of such an operator provide two-sided estimates 
of the eigenvalues of the corresponding restriction of the inverse of the composition of Riesz and 
Bessel potentials, for certain range of parameter a (see, for example, Chen and Song m)- 

Our main goal is the study of the local regularity (modulus of continuity) of the derived weak- 
sense Gaussian solution to equations Sufficient conditions are formulated to obtain the 

mean-quadratic local asymptotic order of the temporal, spatial and spatiotemporal increment 
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random fields, associated with the weak-sense Gaussian solution to equations ©-(H]) (see Theorems 
14.1115.11 and 16.11 below). Specifically, the results derived hold under the condition that the regular 
bounded open domain D is such that the eigenvectors of the Dirichlet negative Laplacian operator 
on D are uniformly bounded. Some examples of domains H, where this condition is satisfied, are 
provided in Section 0 Furthermore, the mean-square Holder continuity in time of the random held 
solution is obtained under some restrictions on the parameter space. While its mean-square Holder 
continuity in space requires the Holder continuity of the eigenvectors of the Dirichlet negative 
Laplacian operator on domain D. The mean-square Holder continuity in space and time directly 
follows, under the above conditions. Also, under such conditions, the sample-path local asymptotic 
orders are obtained immediately from Theorem 3.3.3 of Adler [T] (see Theorem 16.21 below). Note 
that, although the time fractional differentiation in equation © is understood in the weak-sense, 
the Gaussian solution c is Holder continuous in the mean-square sense under the conditions assumed 
in this paper. 

This paper does not adopt the classical framework of diffusion processes characterized by the 
Kolmogorov forward equation or Fokker-Planck equation (see, for example, |26j). In our case, the 
regularized fractional derivative in time, or fractional-in-time derivative in the Gaputo-Djrbashian 
sense, and the Fokker-Planck operator with constant coefficients are applied, in the mean-square 
sense, to a spatiotemporal Gaussian random field for its almost decorrelation in space and time. 
Hence, local self-similarity is observed in the correlation structure in space and in time of the weak- 
sense mean-square Gaussian solution c, as we will prove in this paper. The approach adopted is then 
different from that considered in Chen, Meerschaert and Nane [H], since, in the latter approach, 
the properties of the transition probability densities are investigated, while, in this paper, new 
classes of spatiotemporal Gaussian random fields, displaying local self-similarity, are introduced in 
the weak sense. In particular, the exponents of their local self-similarity are computed in time, 
space and space-time, in the mean-square and sample-path sense. 

Finally, we recall the interest of considering especially spatiotemporal Gaussian random fields on 
Dirichlet regular bounded open domains (see Fuglede [H]), including, as particular cases, bounded 
open C°°- domains, domains with C^-boundary, with Lipschitz continuous boundary, or with 
fractal boundary, among others. Special attention , in the current literature, has been paid to 
the unit ball and the unit sphere, motivated by the analysis of Cosmic Microwave Background 
(CMB) radiation (see, for example, Leonenko and Sakhno [35]; Malyarenko [35] ; Marinucci and 
Peccati m)- In this setting, tensor-valued random fields on the unit sphere are considered for the 
investigation of the combinations Q ± iU, with Q and U respectively representing the linear and 
circular polarization Stokes parameters. 

The outline of the paper is as follows. Preliminary elements and results are presented in 
Section |5| The derivation of a weak-sense mean-square Gaussian solution to equations ©-© is 
established in Section |3| The mean-quadratic local variation exponents in time of the derived 
solution are obtained in Section |3| The mean-quadratic local variation exponents in space are 
given in Section |3| Section [B] then provides the asymptotic local mean quadratic orders in space 
and time. The modulus of continuity of the sample paths of the weak-sense mean-square solution 
to equations ©~© is also derived in this section . Some examples are provided in Section 0 for 
illustration purposes. The extended formulation of the results derived for fractional polynomials 
of the Dirichlet negative Laplacian operator are presented in Section |5| Final comments and some 
open research lines are discussed in Section |3| 

2 Preliminaries 

Some preliminary definitions and results needed in the development of this paper are now intro¬ 
duced. Specifically, some basic results on spectral calculus for self-adjoint operators on a Hilbert 
space are given in Section 12.11 The Mittag-Leffler function is conisdered in Section 12.21 Basic 
elements on fractional Sobolev spaces on a regular bounded open domain are presented in Section 
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2.1 Spectral theory of self-adjoint operators on a separable Hilbert space 

Let us first consider some results on spectral calculus for self-adjoint operators on a Hilbert space. 

Theorem 2.1. (Dautray and Lions, 1990, pp. 119-120 Hd|/ ) Let H be a separable Hilbert space, 

then an injection mapping a exists from the set of spectral families in H into the set of self-adjoint 

operators on H. The following assertions hold: 

Let A be the self-adjoint operator associated with the spectral family ’ where A denotes 

the spectrum of A. The domain of A^ is defined by 

D {A’^) = l^x € H : J X'^’^d{Exx,x) <ooj , k>l. (7) 

For all X G D (A^) , and for all y € H, 

{A'‘x,y)H = j \^d{ExX,y), ( 8 ) 

Ja 

= [ X‘^^d{Exx,x). ( 9 ) 

JA 

If Pk (A) is a polynomial of degree k, then, for all x € D (A^) , and for all y € H, Pk (A) is 
given by 


{Pk{A)x,y)H = [ Pk{X)d{ExX,y). (10) 

JA 

Finally, for a continuous function f on A, the following identities hold for every x € D (f (A)), 
and y G H, 


if i^)x^y)H = f f {X)d{Exx,y). (11) 

JA 

Theorem 2.2. (Dautray and Lions, 1990, p. IjO Let A be a self-adjoint operator in a 
separable Hilbert space H. Lf we denote f the complex conjugate function for f, then D (/ (A)) = 
D{f{A)). Moreover we have {f {A) x,y)H = {x*, f (A) y*), for all x, yGD{f{A)). 

For X G D [f (A)), and y G D {g (A)), then 


{f {A)x,giA)y)H = [ f (X) g (X) d{Exx,y). (12) 

JA 

Furthermore, (/ -|- g) (A) x = f (A) x -\- g (A) x, for all x G D (/ (A)) Gi D {g (A)). 

Finally, if x G D {f (A)), with {g o f) (A) = g {X) f (A), then [g (A) / (A)] x = {g o f) (A) x. 

Theorem IQ is now applied to derive the asymptotic order of the eigenvalues of operator 
(-A)J'" (/—Ad)'’'/^, with, as before, {—A)d representing the Dirichlet negative Laplacian operator 
on regular bounded open domain D. 

Corollary 1. The following asymptotic order holds for the eigenvalues of (/ — Aai)'^/^(—Aa))“/^ : 


A4(-A^)“/2(/-A,,)-r/2) ^ 

fclToo-- =cin,a + j)\D\ P , 

where c(n, a -I- 7 ) is a positive constant depending on n, a and 7 . 


(13) 
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Futhermore, for {(l)k}k>i being the eigenvector system of the Dirichlet negative Laplacian op¬ 
erator {—Ad) on domain D, the following equality holds: 

(-Az5)“/2(/-A,3r/2</.fc = Afc((-Aor/2(/_Az5r/')</)fc, k>l. (14) 

Proof. It is well-known that the eigenvalues { 7 fc(—A£))}fe>i of the Dirichlet negative Laplacian 
operator on domain D C M", arranged in decreasing order of their modulus magnitude satisfy (see, 
for example, Chen and Song my- 

(r (l -I- 

7fe(-AD)-47r - k — ^oo, (15) 


where f{k) ^ g{k) means that limfe_>oo f{k)/g{k) = C, for certain positive constant C. In particular, 
C = 1 in (ITHll . 

From equation (ED in Theorem 12.11 considering f{u) = m“/^(1 -I- we obtain 

Afe ((-Ab)“/2 (/ _ = (7fc(-A,,))“/2 (1 + jki-Ao))^^^. (16) 


Equation m then follows from equations ED and dUD- 

Equation (fHl) is straightforwardly obtained from equation (fTTl) in Theorem 12.11 since in our 
case, i.e., for f{u) = u“/^(l -|- u)'’'/^, for all x,y G H = LF‘{D), with L‘^{D) denoting the space of 
square integrable functions on D, 


f{X)d{E^x,y) = ^(yfc(-A,,))“/^(l + 7fc(-A,,))'^/2 

4'k{u)(l)k {•v)x{u)y{v)dud'v 

I 

= ^ (1 + Jk{-^D)y^^Xkyk: 


fc=l 

x/ 


iDxD 


fc=l 


(17) 


with, as before, {4>k}k>i being the eigenvector system of the Dirichlet negative Laplacian operator 
on D. Specifically, our spectral family is defined in terms of the spectral kernel '/'fc(u)<('fe(v), 

and the spectral measure is given by a point or counting measure with atoms located at the 
eigenvalues. ■ 


2.2 Mittag-Lefiier function 

The weak-sense solution derived in the next section involves the Mittag-Lefiier function. The 
definition of the Mittag-Leffler function, and a two-sided uniform inequality are now considered. 

Definition 1. The Mittag-Leffler function is given by 

= 0</3<l (18) 

(see Erdelyi et al. [12]; Haubold, Mathai and Saxena for a more detailed description of this 
function and its properties). 

Lemma 1. For every j3 S (0,1), the uniform estimate 

l + r(l-^)x - - l + [F(l + ^)]-ia: 

holds over K+ with optimal constants (see Simon m Theorem 4)- 
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2.3 Fractional Sobolev spaces on regular bounded open domains 

The scale of fractional Sobolev spaces is introduced within the spaces S (R”), the space of C°°- 
functions with rapid decay at infinity, and V (R”), the space of C°°-functions with compact support 
contained in R”. The dual of these spaces are respectively the space of tempered distributions, 
S' (R"), and the space of distributions, V (R"). 

For s G R, we denote by H" (R") the space of tempered distributions u such that (l + || u G 

L 2 (R"), A G R". For a regular bounded open domain D in R", we denote 

h" (D) = {u € H''{R") : supp u CD} , (19) 

H" {D) = {f CV {D) : 3F G FT* (R”) such that f = Fd} , (20) 

where denotes the restriction of F to H. With the quotient norm 

iF® [D) is a Hilbert space (see Dautray and Lions, [16], p. 118). 


3 The mean-square Gaussian solution in the weak sense 


The preliminaries given in the previous section are now applied in the derivation of a zero-mean 
Gaussian solution to the stochastic boundary value problem ©-(HI), in the mean-square and weak 
senses. The following result first establishes the suitable range of parameter a and 7 for the 
construction of a Green operator in the trace class, with kernel, the fundamental solution to the 
deterministic problem corresponding to 0 “®- Namely, the following proposition states the ranges 
of parameters a and 7 such that the sequence 


[Ep (-Afc k> 1} 

is in the space of absolute summable sequences, for every f > 0 . 

Proposition 1. For n < a + 

00 

^ Ep (-Afe ((-Ac)“/ 2 (/ - tf^) < 00 , 


for every t > 0. 

Proof. From equation (fT^ . 


hm - ^= c{n,a + -f)\D\ 


( 21 ) 


( 22 ) 


(23) 


Therefore, there exists Fq such that for k > ko, 


(^(-Ai5)“/2(/ - AdT''^^ < 


(24) 


for certain positive constants 0 < Li < L 2 , depending on ko, and a -I- 7 and n. In particular, for 
k > ko, 

1 / I /O', 

1 + [r(l + /3)]-iAfc ((-Ai3)«/2(/ - Aop/^) tP - 1 + [r(l + /3)]-iLiF(“+T')/"t/3’ ^ 

Now, from equations dm) and (E3, applying Lemma [TJ for each fixed t > 0, 
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oo 


fco 

k^l 

oo 

+ ^ Ep(-tPXk({-lXDr/^{I - /XdV/^)) 

k—kQ-\-l 


= M{^,a,j,n) + 

< M{l3,a,j,n) + 

< M{P,a,j,n) + 
= M{P,a,j,n) + 

< M{P,a,j,n) + 


k—kQ-\-l 


oo 

£ Ep[-t^Xk{{-ADr/Hl- 

k—kQ-\-l 

oo 

-fco+1 

^—/3n/{a+~/) ^oo ^ 

-. -/ Ep(u)u'^~^du 

{a + 7)/n Jo ^ 

^-^ra/(a+7) roc ^^^-1 


{a + 'y)/n Jo l + [r(l + /3)] 



du < 00 , 


since 


ko 


M(/3, a, 7 , n) = Ep (-t^Xk ([-AdY^I - ) < oo, 

k^l 

and /p°° i+ [ r^J'+p)]-iu ‘^u < oo, for a + 7 > n. 


(26) 

(27) 


A mean-square Gaussian solution, in the weak sense, to equations ®“0 is formulated in 
Proposition [21 considering D to be a Dirichlet-regular bounded open domain. Note that, in the 
classical theory of boundary value problems, given an open set D with compact closure D in R", the 
classical Dirichlet problem consists of the extension of a given continuous function t/j ■ dD —!■ R to 
a continuous function (j) : D —> R such that (j) is harmonic, that is, satisfies the Laplace equation 
in D. The set D is termed regular if the Dirichlet problem has a (necessarily unique) solution 
for any continuous boundary function tp. For example, every simply connected planar domain is 
regular, but may have a bad boundary, for instance, a fractal boundary (see Arendt and Schleich 
[ 8 ], pp. 54-55; Fuglede [21]). Dirichlet regularity implies that all the eigenfunctions of the Dirichlet 
Laplacian operator on D are bounded continuous functions on this domain that vanish continuously 
on the boundary. This fact will be exploited in the examples given in Section [71 according to the 
conditions required on the eigenfunctions, in the derivation of the main results of this paper. 

In a more general setting, we consider the following definition of Dirichlet-regular bounded 
open domain (see, for example, Brelot [13], p. 137 and Theorem 32, and Fuglede [21], p. 253). 

Definition 2. For a connected bounded open domain D with boundary dD we say that Xq G dD 
is regular if and only if it has a Green kernel such that, for each x G D, 

lim G^(x,y)=0, Vy G D. (28) 

X—>-Xo 


The set D is regular if every point of dD is regular. 






See also Chen et al. m for alternative characterizations of Dirichlet-regular bounded open domains 
in terms of the first exit time in the context of subordinate processes. 

The following result provides a mean-square zero-mean Gaussian solution, in the weak sense, to 
the stochastic pseudodifferential boundary value problem on a Dirichlet-regular bounded 

open domain D. 

Proposition 2. Let c be dehned as 

c(Cx)=/ [ G^(t,x;s,y)e(s,y)dsdy, (29) 

Jo Jd 

where 6 (s,y) is space-time zero-mean Gaussian white noise as given in equation ((21), and 


G^(t,x;s,y) = 

= ^Ep (^-\k {t - (/)fc(x)(/)fc(y), t > s 

k>l 

G^(t,x;s,y) = 0, s > t, (30) 

with, as before, for each k > 1 (see Corollary [IJ 

(-A,3)“/2(/ - A,,)^/2^fe = Afe ((-Az5)“/2(/ _ Aor/^) (31) 

Assume that {(l)k{k>i are uniformly bounded by a constant C{D), depending of the geometrical 
characteristics of the domain D, i.e., 


C{D) = sup (?ifc(x). 

k>l, xGD 


Then, for n < a -I- 7 , c in (ESI) provides a mean-square zero-mean Gaussian solution to problem 
on D, in the weak-sense in the space H°‘ {D). Equivalently, 


ID l 


m.s. 




ilflC {t, x) -h (-A£,)“/^ (/ - AdV^^ c {t, x) 


■!/’(x)dx 


(t, x)'!/;(x)dx, y-ip € ' (D) 


■Q + 7 


'D 


(32) 


where = means the equality in the mean-square sense. In addition, c has covariance kernel given 
by, for all t,s G R+, and x,y G D, 


ptAs p 

i?(t, x;s,y) = E[c(t,x)c(s,y)] = / / G^(t, x; u, z)G^(s, y; u, z)dMdz. (33) 

Jo Jd 

Proof. It is well-known that the solution to the eigenvalue equation 

dP 

—T{t) = -iJiT{t), 0<t<T, (34) 

is given by the Mittag-Leffler function E^{—fj,tP), for any /r > 0, with Ep being introduced in 
equation (ITHl) . Then, for (3 G (0,1), from definition of G^ in equations (l!ini) - (l!TT|) . and the definition 
of the regularized fractional derivative in time ®, 
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Id (^>x;0.y)V'(y)c^y = - ((-Ad)“/2(/- Ad 

xE/ 3 (^-\k ((-Az))“/2(J - Ad)'^/^^ ^fc(x)(/)fc(y)V'(y)dy 

oo 

= - E ((-Az,)“/2(/ - Aor/^) Ep (-Afe ((-A,,)“/2(/ - Ac)^/2) 


X(/)fe(x) / ^fc(y)'0(y)dy 

Jd 

oo 


fc=l 


X(j)k{x)ijjk = -(-Ad)“/^(/ - AnY^'^Qt (V’) 


(35) 


where denotes the integral operator on L'^{D) with kernel (t, x;0,y), for each t > 0. Note 
that, from triangle and Cauchy-Schwarz inequalities, 

°° f)0 

^Ep (-Afe ((-Az,)“/^(/ - AdT/^) t^) M^)^k 
k^l ^ 
oo 

<^Afe ((-A,,)“/ 2 (J-A, 

k^l 

xEp (-Afe ((-A, 3 )“/ 2 (/ - Aor^^) t^) 

OO 

< C{D) Y, \u ({-Aor/^il - AdT^^) 

xEp (-Afe ((-A, 3 )“/ 2 (/ _ Aor^^) i^k 


/c=l 


<C{D)^ 


Y{>^k{{-Anr/Hl-Anp/^)^k? 

\ k=i 


YiEp i-Xk {{-Anr/Hl - AdVG) tP)Y < oo, 

\ fc=i 


where, as before ipk = Jjj (j)k{y)i^{y)dy. Here, 


A Et^'^ ((-Ad)“/ 2 (J - Ad)t'/ 2 ) < oo, 

\ k=i 


since ip G {D), and 


(36) 


Y[Ep (-Afc ((-A,,)“/ 2 (/ - AdV/^) tP)f < oo 

\ fe=i 

from Proposition [TJ since a + 7 > n. 

Applying the regularized fractional derivative in time ([ 6 ]), from equation (13511 . we obtain 
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X = L ro^l L "• •' 


Id r(i — /3) dt Jq 

xe(s, y)-ip(x)dydsdTdx 


J r(i - /3) ^ J ^ ^ J J 
1 


/ ( r(i-/3) / ^ ^ dv^ii;{x)dx 

LIL f(T^ L [s / "■'’®'’ “ - “■ "■ “■ ‘'" 

■* . a/3 


£is,y)dyds 


'0(x)(ix 


^/>(x)dx 

'tp{x)dx 


{-AdT/HI - G^{t,x;s,y) 

Jo Jd 


+ Id I {t,x;s,y)£{s,y)dyds 

= [ ll~^£{t,x)il){x)dx- f 

J D J L 

x£{s,y)dyds\ij){x)dx, (37) 

as we wanted to prove. Here, we have applied that 

OO 

{u,x;u,y) ='^ (l)k{x)(j)k{y) = 6{x-y), Vu > 0, 

fe=i 

with J(x — y) denoting the Dirac Delta distribution on L^{D) such that 

/ iJ(x-y)e(t,y)dy = e(t,x), 

JD 

in the mean-square sense, and in the L^(D)-weak sense. 

Finally, equation (13311 is obtained from straightforward computation of the covariance function 
of c in equation dMl), since for n/2 < a-l-q, G^ defines a Hilbert-Schmidt operator. Consequently, 
its self-convolution defines a covariance operator TZ in the trace class. Thus, its covariance kernel 
R is continuous, and it can be defined pointwise from equation (IMl). ■ 


4 Mean-quadratic local variation in time 

This section provides an upper bound for the mean-quadratic local variation of the temporal 
increments of the mean-square solution c defined in equation (1^^ of Proposition [51 Note that 
although we have showed in Proposition [5] that c satisfies, in the mean-square sense, equation m 
over the test functions in H {D), for a -I- 7 > n, as we prove, in the following result, equation 
(l2^ defines a Holder continuous, in time, spatiotemporal random field c, under a wider range of 
parameter a -I- 7 . Namely, Theorem 14.II below holds for ^ < a -I- 7 , and /3 < 1/2. 

As before, we will consider the sequence of eigenvalues 

Afc((-Az3)“/2(/-Az5)^/2), fc>l, 

arranged in increasing order of their modulus magnitude, with the associated eigenvectors 4>k, 
k > 1 , in the same order. 
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Theorem 4.1. Let c he defined as in of Proposition under the assumption that 

C{D) = supj,>]^ xe-D^fc(^) Then, for f3 < 1/2, and ^ < a + 7 , the following inequality holds: 

E[c{t, x) - c(s, x)]^ < [C{D)]'^g{t - s), (38) 

where 

g{t — s) = O (^{t — , s ^ t, {) < s <t, (39) 

with X Ay denoting the minimum of x and y, for x,y gR. 

Proof. Since Eg{—x) is a monotone decreasing functions with values in the interval [0,1], for 
X £ R+, for 0 < s < t, we obtain 


du 


E[c{t,x.) - c{s,x.)]'^ = f ^/>fc(x) i;, 3 (-Afc ((-Ad)“/^(/-A i>)'^/2^ (t-u)'^) 

•^0 fc=l 

-Ep (-Afe is - uf) 

/ t ^ 

^ cl>liK)E^ (-2A, ((-A^)“/2(/ - An^/^) (t - u)^) 

- fc=i 

n S OO 

< [C{D)f / 5] [Ep (-Afe (i-AnT/^il - AdV/^) (t - uf) 

*^0 1 — 1 


du 




-E^ ((-Afe ((-A,5)“/2(/ _ Any/^) [s - uy) 

nt OO 

+ / E ((-Ad)“/^(/ - Auy/^) it - uf) 

^ S 1 _ 1 




du 


= [cm^ 


nS 

/ ^ [e^ (-Afc ((-A^)“/2(J - Aoy/^) it-u 

+ [e^ ((-Afe ((-Ao)“/2 (/ _ Aor^^) is - u)'^)]' 
2Ep (-Afc (^i-Aor^yi - AdV/^) it - uf) 

xEp ((-Afc (^i-Aor^yi - AdV/^) is - uf) du 

/•t OO 

+ / ^£;/3(-2Afc((-Ac)“/2(/-Az5 

ds 1 _ 1 ^ ^ 

n g CXD 

/ ^ [ep (-Afc [i-Aor^yi - Aoy/^) it - uf) 

Eg ((-Afc ((-A,,)“/2(/ _ AdV/^) is - uf) 

Ep (-Afc (i-AoT/yi - AdV/^) it - uf) 




— uY 1 du 


< [CiDy 
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du 


/ t OO 

^ Ep (-2Afc {t - uf) du 

- k^l 

PS oo 

<[CiD)]^ / ^[i;^(-Afe((-Az,r/2(/-ABr/2)(s-w)/3) 

/O fc=l ^ 

/ i oo 

^ Ep (-2Afc ({-^dT/^I - AdV'^) {t - 

fe=i 

In a similar way to equation (I26|) . from equation (14011 . we obtain 

E[c{t,yd) - c{s,yd)f <[C{D)f [ M{l3,a,-f,n,{s - uf) 

Jo 

+ '^—7- n - / [i^/3 (a;)]^ 2;“+"' 

(a + 7 )/n Jo ^ 


du 


dx 


du 


+ [Cm^ f[Mi/3,a,j,nft-u)f 

J S 

(t — y)-/3"/(“+7) 


(a + 7 )/n Jo 


nOO 

/ i ?/3 (2a;) x°‘+-f~^dx 

Jo 


du 


<[C{D)f [ \M{(],a,j,n,{s-u)f 
Jo 




(s-u)- 


■I 


(a + 7 )/n Jo (1 + [r(l +/3)]-ia;)' 


:dx 


du 


+ [C{D)]'^ f [M{l3,a,-/,n,{t-u)^) 

J S 


{a + f/n Jo 1 + [r(l+ /3)]-i2a; 


dx 


du. 


where 


M(/3,a, 7 , n, {s - u)f = Y.Ep (-2{s - ufXk ((-Az9)“/"(/ - Ad 

k^l 

ko 

M{j3,a,-i,n, {s - u)^) = ^ \Efi {-{s-u)I^Xk {{-Ad)°‘^'^{I - Ad 


fc=l 


Hence, 


i?[c(t,x)-c(s,x)]2 < [C{D)f ki(/3,a,7,n)si-2^ + if2(/3,a,7,^)s'"^"/^“+^^ 


+7^3(/3,a,7,n)(t - s)^ ^ + KflJ,a,^,n)(t - s)'^ /3n/(a+7) 
Thus, when s ^ t, s < t, we have 

E[c{t, x) - c(s, x)]2 < [C{D)]'^g{t - s), 


( 40 ) 


(41) 


(42) 


( 43 ) 
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with 


9{t -s) = 


Ki{P,a,"f,n)s^ + K2{f3,a,"f,n)s^ /3n/(a+7) 
+K3iP, a, 7 , n)(t - a, 7 , n)(t - 

— Q _ g)l-(/3»i/(“+7))A(l-/3)^ ^ 


5 Mean-quadratic local variation in space 

The fractional local asymptotic exponent, in the mean-square sense, of the spatial increments of 
the solution c, derived in Proposition [21 is obtained in the following result. 

Theorem 5.1. Let c be as given in equations for n < a + 'j. Assume that for every 

k>l, 

|,^fc(x-yh)-<^fe(x)| =0(||hr), Ilhll^O, T>0. 

In particular, for each k > 1, and for ||h|| small, 

|(^fc(x -f h) - (^fe(x)l < Cfc||h||'^, 

for certain positive constant Ck- //sup^ Ck = C < 00 , then, for each f > 0, 


| 2 T^ 


Thus, for ||x — y|| sufficiently small. 


yll^O, T>0. 

(44) 


(45) 


where 


g(t) = t^-^^ 


r(i + ; 9 ) 




, t > 0. 


Proof. 

Applying Holder continuity of the eigenvectors, from Lemma [I] we have, for every t > 0, 
A[c(t,x) - c(t,y)]^ 


rt 02 


^ [Ep{-Xk ((-A)“/^(/ - ArJ^) it - uf)M^) 




< 


Ep{-Xk [{-ATJ\i - ArJ^) {t - uf)4>k{y) 

pt 00 

/ ^[</>,(x) - My)? [Epi-Xk {{-ArJ\i - AyJ^) {t - uY) 

•^0 k=l 

pt 00 

C||x - y||2^ / ^ [i?^(-A, ((-A)“/^(/ - A)]/^) {t - uY) 

0 7_1 


du 


nt 00 


<C||x-yf^ / ^ 


r(i + ^) 


0 iYiXk({-A)Y\l-AYY^)vP 


dv 


= + ,.„2T 


tiXk ({-a)Y\i-ayY") 


= Cg(t)\\yi-yY 


(46) 
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as we wanted to prove. Note that, from Corollary[T] for each fixed t > 0, 


for Of + 7 > n. 


9(t) 

/c=l 


r(i + /3) 

A, 


< oo, 


6 Mean quadratic local variation in time and space 

In this section we apply the results derived in Theorems 14.11 and 15.11 to obtain the mean-quadratic 
local variation properties of the spatiotemporal increments of the weak-sense solution c to equations 

®” 0 - 

Theorem 6.1. Under conditions of Theorems and \5.1\. let c be defined in equations i29\) -i31 \} . 
Then, as s ^ t, s,t € (0, T], and ||x — y|| —0, 

E[c{t,x.) - c(s,y)]^ 

< C(iA,T,/ 3 ,a, 7 ,T,n)||(t,x) - (s, 

where 

/A 1/2[(1-^)a(1-/3)A2T]. 

C{D, T, fi, a, 7 , T, n) = 8([C'(il)]2 V Cg{T)) ( - j 

with xy y representing the maximum of x and y, and, as before, x Ay representing the minimum. 
Here, C{D) is introduced in Theorem \4.1[ and C and g{T) are given in Theorem A5.1[ 

Proof. The proof follows from Theorems 14. II and 15. II Specifically, from equations (IHRll - (IMl) . and 
(Bll-dlSI), as s —>■ t, and ||x — y|| —>■ 0 , 

E[c{t, x) - c(s, y)]^ = E[c{t, x) - c(s, x) -f c(s, x) - c(s, y)]^ 

= E[c{t, x) - c(s, x)]^ -f E[c{s, x) - c(s, y)]^ 

+2E [(c(<, x) - c(s, x))(c(s, x) - c(s, y))] 

< E[c{t, x) - c(s, x)]^ -f E[cis, x) - c(s, y)]^ 

-f 2 1 ^; [(c(t, x) - c(s,x))(c(s,x) -c(s,y))]| 


< E[c{t, x) - c(s, x)]^ -f E[c{s, x) - c(s, y)]^ 

-f2 v 'E [{c{t, x) - c(s, x))2 ]^£: [(c(s, x) - c(s, y))^] 

< [C{D)r\t - g|(i -^Wi-/^)+Cg(s)||x-yf^ _ 

+2\/[C{D)]^\t - s|(i-^)^(i-«C'g(s)||x - y||2T 

< [C{D)r\t - + Cg{s)\\K - yf 


+2^ ([CiD)n - + (Cg(s)||x - y|pT) 

< 2{[C{D)]^ V Cg{T)) [it - 

+ ||x - y||2^ + y(|t-,|(i-^Wi-/5 ))V(||x-ypT)2 

i(|t_5|(l-^)A(l-/3)A2T 


<4([C(Zl)]^VCg(T)) 


-|-||x-y||l “+'v 


)a(1-/3)A2T 
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+ yi (|i _ ,|(1 -^)a(1-/3)A2T^2 ^ I'li^ _ y||(l-^)A(l-/3)A2Ty 


( 47 ) 


Under the conditions assumed, 0 < ^1 — A (1 —/3) A2T < 1, hence, we can apply Jensen’s 
inequality for concave function 0 < ^ < 1, obtaining from the last inequality in (I47L 


E[cit, x) - c(s, y)]2 < 4([C(U)]2 V Cg{T)) 

'1\ 1/2[(1-^)a(1-/3)A2T] 

y 


{\t-. 




2 \ (1-^(^)a(1-/3)A2T 


\u 


(|t - s|2 + ||x - y|| 2 )(l-^)^(l-«A 2 T 


1/2[(1-^)a(1-/3)A2T] 

< 4{[C{D)]y Cg{T)) 

X ■(|i_,|2+||x-yf)V2[(l-^)A(l-/3)A2T] 

+ V(|t - S |2 + ||x - y|| 2 )(l-^)^(l-/ 3 )A 2 T 
= 2 C(U,r,/ 3 ,a, 7 ,T,n)||(t,x) - (s, 

as we wanted to prove, with C{D, T, P, a, 7 , T, n) = 2C{D, T, /3, a, 7 , T, n), and 

'i\ 1/2[(1 -^)a(1-/3)A2T] 

C{D, T, P, a, 7 , T, n) = A{[C{D)Y V Cg{T)) ' 


(48) 


6.1 Sample-path properties 

The following result provides the sample path local regularity properties of the mean-square weak- 
sense Gaussian solution c to equations dS])”®- Note that the derived Gaussian solution c is not 
fractional differentiable in time in the strong-sense (see Proposition HI). However, it is Holder 
continuous, in the mean square sense, under the conditions assumed in the previous sections. 

Theorem 6.2. Let c be defined in equations Under the conditions of Theorems \4-l\ and 

15.71 with probability one, the following inequalities hold, as s ^ t, and ||x — y|| 0, 


< 


sup |c(t, x) — c(s,x)p 

|i —s| <5 


^j(1-^)a( 1-/3) iJ^j(l-7^)A(l-/3) 



1/2 


sup |c(t,x) 
l|x-y||<5 
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< + H2S^ 


log| J 


1/2 


sup |c(t,x) - c(s,y)p 

||(t,x)-(s,y)||<5 
< ;^^(1-^)a(1-/3)A2T 

[log Q) 


1 1/2 


(49) 


where Z^Y and X are positive random variables, and Hi, i = 1,2,3, are positive constants that 
could depend on the geometrical characteristics of the domain D considered, like the boundary. 

The proof directly follows from Theorems 14.llfOl and Theorem 3.3.3 in Adler [I], p.57. 


7 Examples 

In the following subsections we consider some special cases of domain D, where the derived results 
can be applied. Specifically, in the examples introduced below, the eigenfunctions of the Dirichlet 
negative Laplacian operator can be explicitly computed, and Theorem 12.11 allows us to define the 
weak-sense mean-square solution to in terms of such eigenfunctions as given in equations 

in Proposition[2] (see, for example, Grebenkov and Nguyen [H]). The conditions required 
in Theorems l4.1l[^?Tll6.1l and l6.21 for the continuity of c in the mean-square sense and in the sample 
path sense, are also verified. 


7.1 Intervals, rectangles, parallelepipeds 

Let us consider the case where D = (0, Li),..., (0, L„) C K", where Li > 0, for z = 1,..., n, 
the method of separation of variables yields the following eigenvectors for the Dirichlet negative 
Laplacian operator: 


; ^n) — si 


sm 


i=l 


7r{ki + l)xt 


(fci,...,fc„) eN", 


(50) 


( k ' “1“ 1 ^ ^ 

with Xki,...,kn = + ■ • ■; j and Xki = — ^ , for z = 1,..., n. In this case, the fundamental 

solution to the fractional space-time pseudodifferential equation 

d^c, 


dtP 


{t,^) = -{-XdT/\I - 


t G R+, X S D — (0, Li),..., (0, L„) 


is given, for a -I- 7 > zz, by 


^ El I / 7r^(fci-|-1)^ 

Lrt —s ( 2:1 , • ■ • , Xn, 2/1 , • ■ • ) Vn) — / , X/p I I y ^ 2 


a/2 




1+ E 


-I-1)^ 


7/2 


{t-sf 


(51) 


n 


sin 


7r(fcj -I- l)xi 


n 


sm 


TT[ki Y l)yi 


, t> s 


and G(t, x; s, y) = 0, s > t. 
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From Holder’s inequality (see, for example, equation (6.4) in Grebenkov and Nguyen [211), 


,...,^^1 (^1 5 ■ • ■ : ^n) I ^ 


fir Ut,L. 

\ 11 * 2’^/2+"7r^/2 ' 

\ i=l 


Therefore, in the previous computations in Theorems 14.11 for /3 < 1/2, we consider 


C{D) = 


2^nl2-\-n^nj2 


Theorem 15.11 also holds, since nr=i ^ (fc^,..., fc„) e N”, are continuously differen¬ 

tiable, and hence. Holder continuous. Theorems 16.11 and 16. 2 1 then follow for (3 < 1/2, and, as before, 
for a -I- 7 > n. 


7.2 Balls 

Let us consider equations @-0 on the ball. That is, 

0/3 „ 

— it,K) = -{-Anr/^{I-Any/^c{t,K), (52) 

teM+, xeT> = {xeM”, ||x|| < i?} = i 3 i{(o), i? > o. 

For a -I- 7 > n, the Green function is defined as 

oo oo h{n,l) 

Gt-s{p,e-,p',e') = Y,J2 E 

l—Q r—1 m—1 

^)0/,r,m(p , ^ ), t ^ S 

and G(t, x; s, y) = 0, s > t, 

where we have considered the spherical coordinates x = {p,0), y = {p',9'), and, as before, Ep is 
the Mittag-Leffler function. Moreover, 


( D4^l,r,m{P: 9^) — j. 0), 

for I G N, with m = 1 , 2,..., h{n, 1), r S N* and h{n, 1) = number of spherical 

harmonics. Here, the eigenvalues are given by fii^r = (%?) ’ eigenvectors (j)i^r,m(,P^9) = 

Jl+n^ (^/+i^,r.fi) Si,m{9) are defined in terms of the Bessel function of the first kind of 
order v, and the orthonormal spherical harmonics on the sphere of radius one, Si^m{9)- Note 
that ci^r,m is the normalizing constant, and is the rth positive root of 

From Holder inequality, since {(f>i^r,m} are normalized in the space of square integrable functions 
over the ball of radius R, i.e.., ||^!>/,r,m ||2 = 1, we obtain 


|</;.r,m(p,0)| < ||</;.r,m||l < Hi ,r,m II2 




(see, for example, equation (6.4) in Grebenkov and Nguyen [24]). Thus, in the previous compu¬ 
tations in Theorem 14.11 we can consider C{D) = |i3ij(0)|^/^. Theorem 15.11 also holds, since Bessel 
functions of the first kind and order //, on a closed interval, and the orthonormal spherical har¬ 
monics on the sphere of radius one, are Holder continuous. Theorems 16.11 and 16.21 then follow for 
/3 < 1 / 2 , and, as before, for a -|- 7 > n. 
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Circular annulus 

Let us now consider, for a + 7 > n, and P < 1/2, 

o/3„ 

x) = -(-Ad)“/^(/ - Any^^cit, x), (53) 

t G R+, D = {{xi,X 2 ) G Rq < \x\ < R}. 

In polar coordinates, xi = rcosip, X 2 = rsini^, the Laplace operator A admits the expression 

, 92 Id 1 92 

A =- 1 - 1 -. 

r dr dip"^ 

The fundamental solution (in polar coordinates) of 


x) = -(-Ad)“/^(/ - Aoy^'^cit, x), 

with 

1_9 1 ^2 \ 

^ \ r i9r dip"^ ) ’ 

is then given by 


00 00 2 


Gt-s{r,(f,r',(p') = XI X! XI(- (^2 fc/i?2)“/^ (1 + {al^f^/R^)y^^ 


n—0k—11—1 


XUn,k,lir, p)Un,k,lir',(p'), 


where 


Un,k,l{f, P) — Jn{TOin,k/R) + Cn^kYnifOin^k/R) X g{Yl(npj 1=2 (71^0) 


with Jn and T„ being the Bessel functions of the first and second kind, and the coefficients an,k 
and Cn,k being set by the boundary conditions at r = i?o, and r = R. 

0 — 7^ ['^ni^n.k') T Cn,kY^(^Cyji,k^] Y h \ jyotn,k') T ^n,kYni^<Xn^kl\ 

K 


0 = - 


^n,k 


i?c 


R 


7/ / '"'J 1 I / “U 

'Jn I ^n^k I Y Cji,k^r) 1 ^n.k r-. 


Jn ( Ctn,; 


Rq 

'“r 


n.kkji 1 OCn^k 


i?C 


Ro 


R 


Y Cji^kYn ( Oiji^k 

H 


Using Holder’s inequality, 


\Un,k,l{r,(p)\ < 



r( 2 ) 


-Rq IliTT, 


,k,l\\2, 


where 


yy" [«, + hV - rf) vipm 


with 


Vn,k{r) = Jn{an,kr/R) Y Cn,kYn{an,kr/R). 


(54) 

(55) 
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From equations o and (EH), the uniform upper bound for the eigenvectors of Dirichlet negative 
Laplacian operator on circular annulus is then given by 


where 


h'^R^ 

< Ml 

Op' . 

n,k 

vl,k{R) 

VI 


since Bessel functions of the first and second kind are uniformly bounded for 0 < Rq < r < R. 
Furthermore, Bessel functions of the first and second kind on a closed bounded interval, as well 
as sine and cosine are also Holder continuous. Thus, Theorem 15.II holds. Summarizing, Theorems 
14.1115.1116.11 and l6.2l hold for a + 7 > n, and P < 1/2. 

Elliptical annulus 

In elliptic coordinates, xi = a cosh r cos 0, X 2 = a sinh r sin 0, the Laplace operator adopts the form: 

a2(sinh^ r + sin^ 9) \dr^ 39^ ) ’ 

where r > 0, and 0 < 0 < 27r, are the radial and angular coordinates, and a > 0 is the prescribed 
distance between the origin and the foci. An ellipse is a curve of constant r = R so that its points 
(xi,X2) satisfy x\/A^ + = 1, where A = a cosh i?, and B = asinhii are the major and 

minor semi-axes, and R denotes the radius of the ellipse. The eccentricity e = aj A = Ij coshi? 
is strictly positive. The interior of an ellipse is characterized by 0 < 0 < 27r and 0 < r < R. 
An elliptical annulus, the interior between two ellipses with the same foci, can be characterized in 
elliptic coordinates (r, 9) with Rq < r < R and 0 < 0 < 27r. 

In elliptic coordinates, the variable separation method, u{r, 9) = g{9)f{r), leads to the following 
equations, after considering that the two differential equations in 9 and r are equal to a constant 

K, 


g"{9) + {K-2qcos29)g{9) =0 (56) 

f"{r) — (k — 2qcosh2r9)f{r) = 0. (57) 

These equations are respectively known as the Mathieu equation and the modified Mathieu equa¬ 
tion, where q = Aa^/4, and the parameter n is called the characteristic value of Mathieu functions, 
whose values lead to a real integer value of the characteristic exponent v of the solution defined 
according to Floquet’s theorem. The two linearly independent periodic solutions of equation ()56|) 
are known as the angular Mathieu functions, and they are respectively denoted as Ken{9,q) and 
sen+i{9, g), n = 0,1, 2,.... That is, we consider k such that the characteristic exponent u satisfies 
iz(k, q) G 1, leading to the referred angular periodic Mathieu functions. There are two linearly inde¬ 
pendent oscillatory radial Mathieu functions of the first kind, solution to equation (1571) . respectively 
denoted as MKn \r, q) and MKn \r, q), corresponding to the same k as Ken{9, q). In addition, there 

is two linearly independent oscillatory radial Mathieu functions of the second kind q) and 

( 2 ) 

Msl^J_^{r,q) corresponding to the same s as sen+i{9,q) (see, for example, Gutierrez-Vega et. al. 
m)- Thus we have four families I = 1,2, 3,4, of eigenfunctions of Laplacian operator in an elliptical 
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domain: 


Unkl{r,0) = Keni0,qnkl)MK^^'>{r,qnkl) 
Unklir,0) = nen{ 0 ,qnk 2 )MK^^\r,qnk 2 ) 
Unki{r,0) = se„+i(6<,g„fc3)Afs^+i(r, gnfcs) 

Unkl{r,0) = sen+i{0,qnk4)Msl^l^{r,qnk4)- 


(58) 


For the elliptical annulus with Dirichlet boundary conditions having radius 0 < Rq < R, there 
are eight individual equations defining the parameter q for each n = 0,l,2 ,..., 

M{R, qnkl) = 0; MkI^\R, qnk2) = 0 

Ms^^l-i^{R, qnks) = 0; Ms^^l^{R,qnk4) =0 

MK^^'^{Ro,qnki) = 0; MK^J^'>{R^,qnk2) = 0 




iRo,qnk3) = 0; Ms^^ji(i?0,gnfe4) = 0. 


The fundamental solution (in elliptic coordinates) of 

d^c. 


dtP 


(t,x) =-(-Az 4)“/^(/-Az4)'^/Mi,x), 


for a + 7 > n, with -A^ = - e) {& + &)’ ^^^en given by 


oo oo 4 

Gt.s{r,0,r',0') = EEE^/3 (- 

n—0k—1 1—1 

X (l + {[Mnkl/a^]))'’^'^ (t - s)^j 

XUn,k.lir,0)Un,k,lir',0'), 

where {un,k,i} are given in equation (1551) . 

From Holder inequality 

\un,k,i{r, 0)\ < M [7r(acoshi?)(asinhi?) — 7r(acoshi?o)(asinhi?o)], 

where M = sup„_^ ^ max(j. g) Unkiir, 0), since angular Mathieu functions, for 0 < 0 < 27r, and radial 
Mathieu functions of first and second kind, for 0 < i?o < r < i?, are uniformly bounded (see, 
for example, Gutierrez-Vega et. al. m)- The uniform upper bound, in Theorems 14.11 for the 
eigenvectors of Dirichlet negative Laplacian operator on D is then given by 

C{D) = M [7r(acoshi?)(asinhi?) — 7r(acoshi?o)(asinhiio)] • 


Finally, since angular Mathieu functions, and radial Mathieu functions of the first and second kind 
on a closed bounded interval are Holder continuous functions, Theorem 15.II also holds. As before. 
Theorems 14.11 [5Tl 16.II and 16.21 follow for a -|- 7 > n, and /3 < 1/2. 


8 Fractional polynomials of the Dirichlet negative Laplacian 
operator on D 


The results formulated in this paper hold under a more general scenario. Specifically, in equations 
©-®, we can replace {—Ad)°‘/‘^{I — A^)^/^ by a fractional elliptic polynomial of the form 






1=0 


{-AdT/\I - AdV'^ 


(59) 
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of degree p, and with constant coefficients c; > 0, Z = 0,... — 1, and Cp > 0. Thus, the following 

reformulation of equation m is considered 

xeA (60) 

with the boundary and initial conditions given in and with 

P ((—A£i)“/^(/— A^))'*'/^) being defined in (1551) . Here, as before, e represents Gaussian space- 
time white noise. The next result provides the extension of the previously established statements, 
for equations ©-®, to equation dsni)- 

Theorem 8.1. The following assertions hold: 

(i) For n < p{a + 7), 

00 

(^-Xk (^P <00, (61) 

k^l 

for every t > 0, where, for each k > 1 

P ((-Ab)“/2(j _ AdV^^) = Xk (P ((-Ao)“/"(/ - Azj)'^/")) (62) 

For n < p{a + "f), the weak-sense solution on F[^^°‘~^'^\d) to Ibdl) . in the mean-square sense, 
with boundary and initial conditions is then given by 

c(i,x)=/' / G^(t,x;s,y)e(s,y)dsdy, (63) 

Jo Jd 

where the integral is understood in the mean-square sense, e{s, y) is space-time zero-mean 
Gaussian white noise as given in equation and, for t > s, 


G^(t,x;s,y) = 

= YEp (-Afe (P (^{-AdY^I - Anr/^)) {t - sf) M^)My), 

k>l 

G^(i,x;s,y) =0, s > t, (64) 

with {(/)fe}fe>i and {Afe (P ((-A_d)“/2(/- Ad)^'/^)) satisfying 
(a) For P < 1/2, and ^ < p{a -\- 7), the following inequality holds: 

E[c{t, x) - c(s, x)]2 < [C{D)]'^g{t - s), (65) 

where C{D) is defined as in Theorem EH and 

g{t — s) = O (^{t — s)(^“p(a+T))^(^“^)j ^ s —> t, 0 < s < t. (66) 

(Hi) For n < p(a 7), assume that the uniform Holder continuity of the eigenvectors of the 
Dirichlet negative Laplacian operator holds, as given in Theorem, 1.5,11 considering ||x — y|| 
sufficiently small, 

P[c(t,x) -c(t,y)]^ < Gg(t)||x-yf'^, (67) 

where C is given in Theorem \5. 11 and 

u) _ ti-/3y'_ + _ t>o 
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(iv) For fj < 1/2, and n < p(a + 7), assume that the uniform Holder continuity of the eigenvectors 
of the Dirichlet negative Laplaeian operator holds, as s ^ t, s,t € (0,T], and ||x — y|| —^ 0, 

£'[c(i,x) - c(s,y)]^ 

< C{D,T,l3,a,'y,p,T,n)\\{t,x) - 


where 


C{D, T, /3, a, 7,p, T, n) = 8([C(i^)]2 V Cg{T)) 


1\1/2[(1-74^)a(1-/3)A2T]. 

V 


(v) For P < 1/2, and n < p{a + 7), under the uniform Holder continuity of the eigenvectors of the 
Dirichlet negative Laplaeian operator, as s ^ t, s,t G (0,T], and ||x — y|| —^ 0, 


< 


sup |c(t, x) — c(s, x)p 

\t — s\<6 


log 


-I 1/2 


sup |c(i,x) - c(i,y)|" 
llx-ylK'S 




log(i 


1/2 


sup |c(i,x) - c(s,y)p 
||(t,x)-(s,y)||<<5 

+H 36 ^^~ (<»+”)p)a(1-/3)A2T 


log| J 


1/2 


( 68 ) 


where Z,Y and X are positive random variables, and Hi, i = 1,2,3, are positive constants that 
could depend on the geometrical characteristics of the domain D considered, like the boundary. 

Proof. (i) As, in Corollary [U we apply equation (HU in Theorem 12.11 considering f{u) = 
P(m“/^( 1 + with P{u) = X)r=o given in equation (l5^ . to obtaining 


Afe (P ((-Azi)“/ 2(/ _ Azi)'^/^)) = ^ Q [(7fc(-Azi))“/" (1 + 7fc(-Azi))'^/^]', 


(69) 


1=0 


where, as before, {jki—^D)}k>i denotes the eigenvalues of the Dirichlet negative Laplaeian oper¬ 
ator arranged in decreasing order of their magnitude. Since 


\2/n 




|P|2/" 

(see, for example, Chen and Song m), we obtain 

Afe (P((-Azi)“/2(/_a,,)t'/2)) 


lim 

k — ycx) 


J^{a+-i)p/n 


= d(n,a + 7,p)|P|-P(^+“)/", 


(70) 


(71) 


where c(n, a + 7) is a positive constant depending on n, a, 7 and p. 

Futhermore, equation (HU in Theorem HU also implies the following equality: For fc > 1, 

P ([-XdT/^I - c^k = \k (p ((-A,3)“/2(/ - A,,)'^/2^) (72) 
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for the eigenvector system {4>k}k>i of the Dirichlet negative Laplacian operator (—Ad) on domain 

D. 

From equation (ED, there exists ko such that for k > ko, 

l^kPio^+l)/n < (73) 

for certain positive constants 0 < Li < L 2 , depending on ko, and p{a + 7 ) and n. In particular, 
for k > ko, 

1 

1 + [r(l + /3)]-iAfe {P ((-Ad)“/2(/ - Ad)^/2)) 

< --. (74) 

1 + [r(l +/3)]-iLifcp(“+7)/"t/3 

From Lemma [U in a similar way to Proposition [TJ for each fixed t > 0, 


00 

^ Ep (-t^Xk (P ((-Ad)“/2(/ - Ad)^/^))) 


< M{l3,a,-f,p,n) + 


f-Pnjp{a+'^) 

p{a + 7 )/n 



y^p{a+’y) 


l + [F(l + /3)]-iu 


du < 00 , 


(75) 


Since 


ko 


M{/3,a,j,p,n) = ^Ep (-t^Xu (p ((-Ad)“/'(/ - Ad 

k^l 

fiT i+[r(i+l\]~^ n^'^ < p(^ + 7 ) > n. 

Applying triangle and Caue 
equation (1751) holds, we obtain 


< 00 , 


(76) 


p ( Oi "y ) 

Applying triangle and Cauchy-Schwarz inequalities, since 'll) G H (^); and, for p(a+ 7 ) > n 


^Ep (-Xk (p ((-Ad)“/2(/ - Ad)'^/^)) t^) 

k^l ^ 

CO 

< ^ Afe (P ((-Ad)“/^(/ - Ad)"^/^)) 

k^l 

xEp (-Afe (P ((-Ad)“/"(/ - Ad)^/")) IM^Mk 

CO 

< C{D) Y, Xu (p ((-Ad)“/'(/ - Ad)^/')) 

k^l 

xEp (-Afe (P ([-AdT/^I - Ad)^/')) t^) ■>Pk < 00 , 

(77) 


where, as before, 'i/’fe = Jo 4>k{y)'ip{y)dy. 

Applying the regularized fractional derivative in time ([HI), we then obtain, in a similar way to 
Proposition [U 
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J ^c(t,x)'!/'(x)dx = y ll '^e(t,x)V'(x)dx 

-J J J {t,x;s,y)e{s,y)dyd. 


'(/i(x)dx, 


( 78 ) 


J) ^ Q/ _|_ 'Y ^ 

for every ip a H {D), as we wanted to prove, 
(ii) In a similar way to Theorem 14.11 


£;[c(t,x) - c(s,x)]2 < [(7(1))]^ / M(l3,a,j,p,n,{s -u)^) 

Jo 

(s _ y)-/3"/p(a+7) /■“ 2,p(c+7)“l 


+ - 


p(a + 7)/n Jo (1 + [r(l+ ^)]-ix)' 


dx 


du 


+ [C{D)]‘^ f [M(/3,a,7,p,n, (t-w)^) 
J S 

(^i — y)-/3’T’/p(«+7) r°° 


+- 


■L 


j^p{a + ~f) 


-1 


p(a + 7 )/n Jo 1 + [r(l +/3)]-i2x 


dx 


du. 


where 

ka 

M(/3,a,7,P,n, (s - uf) = Y.Ep (-2{s - ufXk (p ((-Az5)“/2(J - )) 

ko 

M{/3,a,j,p,n, (s - m)^) = ^ (^-(s - ufXk (p ((-Ac)“/2(/ - 

fc=l 


Hence, 

P[c(t,x) - c(s,x)]^ < [^(P)]^ [Ki{l3,a,j,p,n)s~^~'^^ 

+A:2(;9,a,7,p,n)s^“'^"/(P(“+'^^) + Ko{l3,a,^,p,n){t - 
+K4{I3, a, 7 ,p, n)(t - s)i-/3V((a+7)p) , 

Thus, we have 

E[c{t, x) - c(s, x)]2 < [C'(P)] 25 (t - s), 

with, for s —)• t, s < t, 

g(t - s) = O ((t - s)1-(/3"/p(«+7))a(i-/3)'\ ^ 


(79) 


(80) 


(81) 


(iii) Applying Holder continuity of the eigenvectors, from Lemma [1] in a similar way to Theorem 
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[SH for every f > 0 , 


E[c{t, x) - c(f, y)]^ < Cllx - y f 


E 

fc=i 


= C||x-y 


|2T 




_r(i + /3)_ 

Afc (p({-^rj\i-/^v^^))uP 

r(i + /3) 


du 




= C'sWIIx-yf’^, 

as we wanted to prove. Here, for each fixed t > 0, 




gif) = t^ ^E■ 


r(i + /3) 


( 82 ) 

(83) 


ti Afc (p ((-A)“/2(J- A)^/2)) 

for p(a + 7 ) > n. 

(iv) In a similar way to Theorem 16.11 since under the conditions assumed, 0 < ^1 — p(^a+'y ) ) A (1 — 
/3) A 2T < 1, applying Jensen’s inequality we obtain. 




E[cit, x) - c(s, y)]2 < 4([C(P)]2 V Cg{T)) J 

+ \/(|t - s |2 + ||x - y|| 2 )(l-i;Tlf^)^(l-/ 5 )A 2 T 

= 2 C(P,r,/ 3 ,a, 7 ,P,T,n)||(f,x) - (s,y)||('-^)^('-^)^''', 
where C{D, T, /3, a, 7 ,p, T, n) = 2C{D, T, P, a, 7 ,p, T, n), and 

/1\1/2[(1-5;^)a(1-/3)A2T] 

C(D,T,/],a,j,p,r,n)=4(lC(D)]^VCg(T)) f-j 

(v) The sample-path regularity properties follow straightforward from (ii)-(iv), by applying The¬ 
orem 3.3.3 in Adler [I], p.57. 


(84) 


9 Final comments 

Under the conditions assumed in Proposition [TJ a mean-square solution to equations is 

derived in Proposition [2j in the weak-sense on the space (respectively on the space 

H {^): ill general case considered in Theorem IS.lD . In particular, from the results derived, 
we can define a 7J“(“+''')(P)—valued stochastic process {Ct, t € R'*'}, on the basic probability space 
(H,yl,P), satisfying equation ([3]) a.s., i.e., 

(Ct{; a;), E + (/ - 

= (7y^£(-,cc),V’(-))^,^^^, a.s, V^e77“+^(P), teR+. (85) 
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(Note that similar assertions hold for the derived solution to equation in Theorem 18.11) . In 
this derivation, the orthogonality in of the random components of e is applied, i.e., 

we have applied that E[e{t,x)] = 0, and iil[e(t,x)e(s, y)] = S{t — s)S(x — y), for t,s G IR+, 
and x,y G _D C K". It is well-known that this property holds for any white noise measure on 
L^(]R_i_ X D), beyond the Gaussian case. Furthermore, the fractional integration in the definition of 
the driven process ll~^s in equation m is understood in the mean-square sense on a suitable space 
of test functions, as given in Proposition [2l Thus, we have only considered the properties of the 
second-order moments of the distribution of the driven process in equation (|31). Hence, Proposition 
[5] also holds when the Gaussian space-time white noise on L^(IR+ x D) is replaced by an arbitrary 
white noise random measure e on x D). In particular. Levy noise can be considered. In 

that case, Theorems 14.11 [^?T] and IHTTl respectively provide the Holder continuity, in the mean-square 
sense (i.e., the continuity of the second-order moments), in time, space, and space and time of the 
weak-sense solution, defined by integration with respect to Levy noise measure drj, as 

cit,x)=f [ G^{t,x]s,y)dr]{s,y), (86) 

Jo Jd 

with, as before, being defined in (|^ . for the case of Proposition [2j and in (l64ll . for the 
case of Theorem O In both cases, we can interpret the integral (1861) as a multiparameter ltd 
integral with respect to n -I- I-parameter Levy process r], since, for a -I- 7 > n, and for each 
t G R+, Gf* defines a trace operator Qt on L?{D), as proved in Proposition [T] (see, for example, 
Lpkka, 0ksendal and Proske [53] j for an alternative interpretation and derivation of solutions 
in that Levy noise case, in terms of functions with values in the Kondratiev space of stochastic 
distributions). Summarizing, the derived results provide the characterization of the second-order 
regularity properties of the weak-sense solution to equations (I3|)-(|3|) (respectively, to equation (15(1)) 
in Theorem IQ) . For the non-Gaussian case, further research should be developed in order to 
obtain the distributional characteristics of (I86L beyond the second-order moments. This subject 
will be considered in a subsequent paper. Note also that Theorem 16.21 on the characterization 
of the sample-path regularity properties of the weak-sense solution to equations (I3])-(|3|), in the 
mean-square sense (respectively, (v) of Theorem 18.11) only holds for the Gaussian case. 

The authors have recently got new results on fractional-in-time and multifractional-in-space 
stochastic partial differential equations. Such results also appear in the Journal Fractional Cal¬ 
culus &: Applied Analysis, Vol. 19, pp. 1434 1459 , DOI: 10.1515/fca-2016-0074, and 
is available online at http://www.degruyter.eom/view/j/fca 
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